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Abstract 

We study exact stationary properties of the one-dimensional Kardar-Parisi-Zhang 
(KPZ) equation by using the replica approach. The stationary state for the KPZ 
equation is realized by setting the initial condition the two-sided Brownian motion 
(BM) with respect to the space variable. Developing techniques for dealing with this 
initial condition in the replica analysis, we elucidate some exact nature of the height 
fluctuation for the KPZ equation. In particular, we obtain an explicit representation 
of the probability distribution of the height in terms of the Fredholm determinants. 
Furthermore from this expression, we also get the exact expression of the space-time 
two-point correlation function. 

1 Introduction 

Surface growth phenomena have attracted much attention in both science and technology, 
which involve fundamental understanding of the roughness of the surface and applications 
to the design of new materials. Being typical nonlinear and nonequilibrium phenomena, 
they have been one of the main topics in nonequilibrium statistical physics. In 1986, 
Kardar, Parisi and Zhang proposed a prototypical equation which describes surface growth 
with local interaction [1]. Its one-dimensional version reads 

dh(x, t) _ X fdh(x, t)Y + Mzt) + ^ (a . ty (L1) 



dt 2 \ dx J dx 

Here h(x,t) represents the height of the surface at position x G M and time t > 0. The 
first term represents the effect of nonlinearity and the second one describes the smoothing 
effect of the surface. rj(x,t) represents randomness described by the Gaussian white noise 
with covariance, 

(t}(x, t)rj(x', £')> = 8(x - x')5(t - t'). (1.2) 
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The parameters A, u, D determine their strengths. 

An important feature of the KPZ equation is that it can describe self-affine properties of 
the height fluctuation. In one- dimensional case, it has been shown by the renormalization 
techniques that the height fluctuation grows like t 1 ^ 3 as t goes to infinity [1]. This growth 
exponent 1/3 appear also in many models of surface growth and characterizes the KPZ 
universality class [2, 3]. 

Our understanding of the KPZ universality class has progressed by finding mathemat- 
ical connection between the probability distribution of the height and the random matrix 
theory. The first breakthrough was brought not in the KPZ equation but in solvable dis- 
crete models belonging to the KPZ universality class such as the totally asymmetric simple 
exclusion process (TASEP) and the polynuclear growth (PNG) model. In [4], the author 
obtained the exact solution of the current distribution for the TASEP with the step initial 
configuration, which corresponds to the height distribution for the curved growth in the 
language of growth process. It turned out that it converges in the long time limit to the 
largest eigenvalue distribution of the Gaussian Unitary Ensemble (GUE) in the random 
matrix theory. This distribution is called the GUE Tracy- Widom distribution [5]. The 
result in [4] has been generalized to the case of the other initial configurations [6, 7] and 
it has been found that the current distributions depend on the initial conditions even if 
they have the common scaling exponents 1/3. For example, in the case of the alternating 
initial condition for TASEP, which corresponds to the flat initial configuration in a growth 
process, the limiting current distribution converges to the Gaussian Orthogonal Ensemble 
(GOE) Tracy- Widom distribution, which is the largest eigenvalue distribution of the other 
random matrix ensemble called the GOE [6, 8, 9]. Since then, the studies on the height 
distributions have been the subject of active investigation [10, 11, 12]. 

During the last few years, the studies on the KPZ equation and KPZ universality have 
entered a new stage [13, 14]. First, an experiment using the liquid crystal turbulence were 
performed in [15, 16, 17]. The authors succeeded in obtaining not only the exponent 1/3 
but also the GUE/GOE Tracy- Widom distributions as well as their finite-time corrections. 
Second, there has been much progress on the theoretical study: An exact height distribution 
for the KPZ equation has been found in [18, 19, 20, 21]. This first exact solution was 
obtained for the narrow wedge initial condition, from which the surface grows in parabolic 
shape. The analysis was based on the result of the current distribution of the ASEP 
in [22, 23] and the fact that the KPZ equation is obtained as a weakly asymmetric limit 
of the ASEP [24]. Afterwards the same technique was applied to the case of the half 
Brownian motion initial condition [25]. A remarkable feature of this solution is that this 
distribution describes the universal crossover between the Edwards- Wilkinson [26] and the 
KPZ universality classes. As t — > 0, the growth exponent becomes 1/4 and the distribution 
converges to the Gaussian, which implies that the nonlinearity in the KPZ equation does 
not become effective and thus the system belongs to the Edwards- Wilkinson universality 
class. As t — > oo, on the other hand, the exponent becomes 1/3 and the distribution 
converges to the GUE Tracy- Widom distribution. Furthermore it turned out that this 
crossover is universal for weakly driven growth [27, 28]. 

As explained above, these achievements in the KPZ equation provide us not only exact 
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information but also the important physical findings. Thus it is of great importance to 
develop the exact approach in such a way that it can be applied to other interesting 
situations. In particular, a most challenging problem would be an application of the exact 
solution to the analysis on the space-time correlation in the case of the stationary situation, 
which is the main subject of this paper. For this purpose, the replica method is suitable. 
The replica method has been established as a powerful tool for disordered systems in 
statistical mechanics. The application to the KPZ equation was first proposed in [29] and 
refined in [28, 30, 31]. It utilizes a connection between the iVth moment of the exponential 
height (N replica partition function) and the dynamics of one-dimensional N particle Bose 
gas system with attractive interaction which can be solved by the Bethe ansatz [32, 33]. 
Using the exact knowledge to the N replica partition function obtained by the Bethe ansatz 
technique, we analyze the generating function of the N replica partition function. Although 
it is in fact a divergent sum, a resummation technique developed in [28, 30, 31] allows us 
to find Fredholm determinant representation of the generating function from which the 
exact height distribution is readily obtained. The method is so powerful that not merely 
the same expression was rederived for the narrow wedge initial condition [28, 30, 31], but 
one can generalize it to the case of other initial conditions, such as the flat [34, 35] and 
half-Brownian motion [36] initial condition. The multi-point distributions have also been 
studied [37, 38, 36] though it includes a decoupling assumption. Furthermore the idea of 
finding the Fredholm determinant structure of moment generating function has become 
the basis of the study on the Macdonald process [39] , which is a larger class of integrable 
stochastic process including the KPZ equation, a directed polymer model [40], etc. and 
deals with the rigorous version of the replica analysis. For more recent developments, 
see [41, 42]. 

In this paper, we consider the stationary situation of the KPZ equation by the replica 
method. It is known that the stationary state is given by the Brownian motion (BM) [43] 
with respect to the position (x) axis i.e., the stationary state is realized when we prepare 
the initial height difference h(x, 0) — h(0, 0) as the two-sided Brownian motion. Since the 
system is translational invariant, one can set without loss of generality the initial height a 
the origin to be zero, i.e., h(0,0) = 0. Then the initial condition considered in this paper 
is given by 



where a = (2z/)~ 3//2 AD 1//2 and B±{x) are two independent standard BMs with expectation 
E[B(x)} = and covariance Cov[B(x)B(y)] = min{x, y}. For this initial condition we 
consider the height distribution function 



and obtain its explicit expression in terms of the Fredholm determinant. From this expres- 
sion, we also obtain the stationary two point correlation function, 




x < 0, 
x > 0, 



(1.3) 



Prob (h(x, t) < s) 



(1.4) 



((h(x,t)-{h(x,t))) 2 ), 



(1.5) 
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which is one of the most fundamental quantities which detects the space-time correlation 
in statistical mechanical systems. Note that (1.5) describes the correlation between h(x,t) 
and h(0, 0) although h(0, 0) does not appear in this expression since we set h(0, 0) = 
as explained above. We also show that in the long-time limit, scaled forms of these 
functions (1.4) and (1.5) converge to the ones obtained in [7] and [44, 45, 46, 47] respectively. 
A part of these result in this paper was already reported in [48]. The purpose of this paper 
is to provide complete descriptions of our results including their derivations. 

The paper is organized as follows. In Sec. 2, we will give precise definitions of the 
quantities discussed in this paper and our results. In Sec. 3, the details of the replica 
analysis will be explained. In Sec. 4, we will obtain the exact solution for the generating 
function of the replica partition function and the height distribution function for the two- 
sided Brownian motion initial condition with drifts. By taking the zero drift limit, our 
final goal, the height distribution for the stationary state, is accomplished, which will be 
explained in Sec. 5. The long-time limit will be taken in Sec. 6. The last section will be 
devoted to the conclusion. 

2 Model and main results 
2.1 Height distribution function 

2.1.1 Two-sided Brownian motion initial condition with drifts 

One of the goals in this paper is to obtain an exact expression for the distribution of the 
height (1.4) with the two-sided Brownian motion initial condition (1.3). For this purpose, 
it is convenient as a first step to consider a slightly more generalized initial condition where 
we add the drift terms to (1.3), 



where v± indicate the strengths of the drifts. Fig. 1 depicts this generalized initial condi- 
tion (2.1). At first we set the drifts to be positive so that one can take an average over the 
Brownian motion initial condition (see (3.15) below). Afterwards, the condition for the 
drifts will be relaxed gradually as our analysis goes on (see (3.35) and (4.21)). In Sec. 5, 
we finally take the limit v ± — > which corresponds to the stationary situation. 

The macroscopic shape of the height h(x,t) for large t is obtained by solving (1.1) 
where only the nonlinear term is taken into account and other terms (the diffusion and 
noise terms) are ignored. For the initial condition (2.1), it becomes 




(2.1) 




(2.2) 
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h(x,0) 




Figure 1: Two-sided Brownian motion initial condition with drifts v± 



where a is given below (1.3). We are interested in the fluctuation property near the 
origin x = around this macroscopic profile. From the KPZ scaling, one expects that the 
fluctuation of the height scales like Oit 1 ^) and nontrivial correlations are seen in the x 
direction with scale 0(t 2 ^ 3 ) when t is large. Let us define a parameter j t which scales as 
0(t^ 3 ) and a rescaled space coordinate X by 

7t= {a 4 ut)\ s = ^£* (2.3) 
We can introduce the scaled height H t (X) by 

± h (M±A = _lt_ JtX i + JtHt{X) . (2 .4) 



2v V a 2 J 12 

Here the term — / -f t X 2 in RHS comes from the parabolic profile (the middle equation) 
in (2.2). Let us define the probability distribution of H t (X), 

F v+>v _, t {s;X) = Prob (J-h ft^,*) + \ + ^X 2 < 7t sJ = Pmb{H t {X) < s). (2.5) 

In order to get an exact expression for this quantity, we now introduce the generalized 
height h(x,t) defined by 

h(x,t) = h{x,t) + X , (2.6) 

where the initial height of the overall surface x is set to be a random variable such that e x 
obeys the inverse gamma distribution with parameter v + + i.e. the probability density 
function of e x is given by 

= tv !l. ^ ^ H ^W ( 2 - 7 ) 

Although this definition (2.6) may look artificial, we will find that it has more tractable 
mathematical structure compared with h(x,t) itself (see Proposition 1 below and Sec. 4.2). 
We also define the probability distribution F v+iV _ tt (s;X) in a similar way to (2.5) 

F v+ , v _, t (s;X) = Prob f-^h (j^,t) + y| + r VtX 2 < its^j = Prob(H t (X) + ^L< S ). 

' (2-8) 
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Our result for F v+>v _ tt (s; X) and F v+tV _ jt (s; X) are summarized as follows. 
Proposition 1 

^v-A*) = r , rr^ v,^). (2-9) 

r(f + + f_+7 t a/as) 

/oo 
due- elt(s ~ u) g lt (u]X). (2.10) 
-oo 

Here g lt (u;X) is expressed as a difference of two Fredholm determinants, 

<? 7 >; X) = det (l - P U (B^ - P*.)P U ) - det (l - P U B^P U ) , (2.11) 
where P u is a projection operator on (u, oo), and 

POO 

<(£i,6) = / dyA? r 
Jo 

poo 

+ / dy 



e ity _ 1 



1 X X\ A . r / 1 X X 

7t 7t 7t/ V 7t 7* 7t 

1 X X\ A . r / 1 X X 

Ai r 6 + — ,«+ — ,«_ + —] Ai r & + — ,v- + —,v + — 

It It It J \ It It It 



1 X X\ A . r / 1 X X 
-Ai r Ki - y, — , u+ , v- + — Ai r £>-?/, — , u- + — , v+ 



It It It J \ It It It, 

(2.12) 

= Ai? ^, --,«- + Aif + - - Y (2.13) 

\ 7i It It J \ It It It J 

and the deformed Airy function Aip(a, 6, c, d) is defined by 



Here T Zp represents the contour from —oo to oo and, along the way, passing below the pole 
z p = id/b. 

We obtained this result by the replica method. The details of this approach and 
the derivation of Proposition 1 will be explained in Sees. 3 and 4. In (2.9), the factor 
— — r ^ v++v l), , , comes from the Laplace transform of the pdf of e x and is defined by 

T(v + +v-+-y t 1 d/ds) J 

the Taylor expansion. This equation represents the relation between F v+!V _ it (s;X) and 
Fv + ,v-,t(s~, X), which will be discussed particularly in Sec. 4.2. 

Note that the distribution function (2.10) is expressed as the convolution of the Gumbel 
distribution and the Fredholm determinants. Such a mathematical structure is common to 
the cases of the narrow wedge initial condition i.e. h(x, 0) = —\x\/5, 5 — >■ [18, 19, 20, 21] 
and the half Brownian motion initial condition, i.e. h(x, 0) is BM for x > while the 
narrow wedge for x < [25, 36]. (This is the reason why we introduce the generalized 
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height (2.6).) Moreover, (2.10) can be regarded as a generalization of the probability 
distributions for the narrow wedge and half-BM initial conditions: Noticing the relations 

(1 X X\ ( 1 X x\ 
x,-,v+ , v. + — AiM y, -, v_ + — , v + = M(x)M(y), (2.15) 
7* 7t 7t/ V 7t 7t 7t/ 



X X\ . . r ( 1 X X 



lim Ai r [x, — ,v+ ,v- + — Ai r y,—,v- + — ,v+ 

\ 7i It It J \ It It It 



= Ai r (x,f ! -^Vi r ^ 1 >--h (2-16) 
where 



It It J V It It 



1 f°° ,3 1 

Ai r (a,6,c) = — / dze faB+i T (2.17) 
27r y_ oc r (-z6z + c) 

Ai r (a, 6, c) = — / ^e i2a+i ^r (ibz + c), (2.18) 

we easily find the function (2.10) becomes that for the narrow wedge case as both v + and 
V- go to infinity and that for the half Brownian motion case as V- (v + ) goes to infinity (0). 



2.1.2 Stationary limit 

Here we consider the stationary limit where both drifts v± in (2.1) go to 0. When we take 
this limit for the expression (2.9), we notice that the Gamma function factor T(v + + 1>_) 
in (2.9) becomes divergent while the factor F v+;V _ ;t (s) vanishes since the term x i n (2-6) 
diverges. Thus we have to analyze carefully the limiting behavior of each factor, which will 
be discussed in Sec. 5. We obtain the following result: Let us introduce the scaled drifts 
u± = v±l r ) t and consider the quantity F Wtt (s) defined by 

F w>t (s;X)= lim F u+/yt ^_ /ytit (s;X). (2.19) 

UJ— — > — CO+ 
<jJ-\-=W 

Note that the situation u + + uj_ = corresponds to the case of flat macroscopic profile 
(see Fig. 1 and (2.2)) and the parameter w determines the slope of the surface. The case 
w = corresponds to the initial condition (1.3). In terms of the functions 

/ d\e aX A&( X + U + X,-,1 + -,1 + - ] , (2.20) 
Jo V It It It) 

(2.21) 

F Wjt (s;X) is written as follows. 
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Theorem 2 



F w ,t(s; X) = — d/ g f° due-^ s ~ u) U t (u; X) - X)) . (2.22) 

r(l + 7 i d/ds) J^oo V / 



Here v Wyt {u;X) and u^ t (u;X) are given by 

v w ,t(u; X) = det (1 - A w - X - w +x) L w - X ,-w+x(u) + det (1 - A w _ x ,-w+x - D w - X ,-w+x) , 

(2.23) 

i/g(u; X) = det (l - A^_ w+X ) L^_ w+X (u) + det (l - A^_ w+X - D^_ w+X ) , 

(2.24) 

respectively. The kernels A a b , L a ^ and the function D a b are represented as 

= C t (6) r dyM T T U + y, 1 + 1 + -) Ai? f 6 + y, -, 1 + -, 1 + - Y (2.25) 
A V 7t 7t 7t/ V It It It) 

L a ,b(u) = ~r / dxC t (x)B aAu (x)B b:ajU (x), (2.26) 

a + o J -oo 

Am>(6, 6) = (4,,&C t £ 0iM ) (£l)£ 6l «,„(6) (2.27) 

where the functions A^\, and D^\(u) are also defined in the same way as (2.25), (2.26) 
and (2.27) respectively with C t (x) (2.21) replaced by Cf\x) := J? t T_ x — 8(x). 

Although in the limit b — > —a, the first term in (2.26) becomes divergent, we find the 
second term compensates for the divergence and obtain the following expression for (2.26), 

L a ,- a (u) = + u - a 2 - 1 + f°° dxC t {x) (B { X,u( x ) B S-a,u(x) + B^_ a>u (x)B^ u (x)) 

It J-OO V 7 

/oo 
dxC t (x)B^_ a!U (x)B^ u (x), (2.28) 
-oo 

where 7 is Euler's constant and 

p a 3 /3-(x+u)a 

A 1 -) 



r(^ + i)' 

/ rfAe aA AiM x + m + A,-,1 + -,1 + - ] (2.29) 
Jo V 7t 7t It) 



are the first and second term in (2.20). 

The expression in Theorem 2 is suitable for numerical evaluations of the distribution 
function. By contrast Proposition 1 is not because we need precise estimations for both the 
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divergence of the Gamma function factor and the convergence of the factor F v+;V _ !t (s; X) to 
zero in the stationary limit. Fig. 2 shows the pictures of the probability density function 
for w = X = 0. These are obtained by discretizing simply the Fredholm determinants 
in (2.23) and (2.24), which is sufficient for our purpose of getting those pictures. A more 
elaborate approximation technique for a Fredholm determinant has been developed recently 
in [49, 50, 51]. For accurate estimations of statistical quantities of our height distributions, 
this technique would be useful. 

The distribution function obtained in this paper appears not merely in the KPZ equa- 
tion but also in some time regime in various stochastic models [24, 27, 28]. One typical 
example is the asymmetric simple exclusion process (ASEP). The ASEP is a stochastic 
many-particle system where each particle hops to the right (left) neighboring site with rate 
q (p) but it cannot hop to the occupied sites, which leads to the exclusion effect. Now we 
consider the ASEP on the one dimensional lattice Z and assume that at t — 0, particles 
obey the two-sided Bernoulli initial condition with density 1/2, where a particle occupy 
each site with probability 1/2. Note that such a Bernoulli product measure is a stationary 
state of the ASEP. Let N(x, t) be the (integrated) particle current i.e. difference between 
the numbers of the particles which pass from the site x to the site x + 1 and from the site 
x + 1 to the site x up to time t. 

In [24], it was shown that 

lim {2N(e~ 1 x, e" 2 t) + te' 3/2 /2) = hi^(x, t) + ^ (2.30) 

where e := q — p is the difference between the right (q) and the left (p) hopping rates. 
The function hi ^x, t) is the (Cole-Hopf) solution to the KPZ equation with parameters 
v = 1/2, A = 1, D = 1 and with the initial condition given in (2.1) with v + = V- = 0. 
This relation indicates that the height of the KPZ equation well approximates the current 
in the ASEP with the asymmetry q — p being small, which is called the weakly ASEP 
(WASEP). Moreover it has been known recently that the KPZ equation also describes the 
dynamics of many other weakly driven growth processes in some time regime [27, 28]. The 
dots in Fig. 2 represent the Monte Carlo simulation for the particle current fluctuation 
yfl (2N(0, e- 2 t) + te~ 3/2 /2) for the WASEP with the two-sided Bernoulli initial condition 
mentioned above. 

2.1.3 Long-time limit 

We can also obtain the long-time limit of the height distribution function. We define 
F w (s; X) as lim^oo F w>t (s; X) and obtain the following result. 

Theorem 3 

F w {s-X) = ^-v w {s-X). (2.31) 
where u w (s; X) := lim^oo v Wjt (s; X) has the following expression. 

u w (s; X) = det (1 - A) £„_x,-„+x(s) + det (1 - A - V w _ x ,_ w+X ) . (2.32) 
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Figure 2: The probability density function of the height in the stationary situation 
dF w=0)t (s; X = 0)/ds for j t = 1,2,3. The dashed curve is dF (s; 0)/ds. The dots (■) 
indicate the Monte Carlo simulation for current distribution for the WASEP. We set the 
right hopping rate q = 0.575 and the left rate p = 0.425. (i.e. e = q — p = 0.15), 3951 
Monte Carlo steps (corresponding to j t = 1), and 10000 samples. 



Here A, P 0i & and C a ,b are the long-time limit of A a ^, D a ^ and L a ^ defined in Theorem 2 
and are given by 



/•OO 

A{x, y) = P {x) / d£Ai(x + £)Ai(y + f), 

J u 



£a,b(s) 



1 — a — b 



dxB a:S (x)B b!S (x), 



a + b j Q 

V a , b {x,y) = (AP B a , s ) (x)B bta (y), 

3 r°° 
B ajS (x) = eT-^ 8 '" - / d\e aX M(x + s + A). 

Jo 



(2.33) 

(2.34) 
(2.35) 
(2.36) 



Note that det (1 — A) in (2.32) is nothing but the GUE Tracy- Widom distribution. The 
expression for C a _ corresponding to (2.28) is given by 

dx (b { X(x)BW(x) + B^(x)B%(x)) - J dxB$(x)B { X(a 

(2.37) 



where B^l(x) (B^l(x)) represents the first (second) term in (2.36), i.e. 



£«(*) = eV-(^)* B?]( 



a,s 



dXe aX A\(x + s + X). 



(2.38) 
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This function F w (s) has already appeared as a limiting distribution for other stochastic 
processes such as the PNG model [7] and the TASEP [46] in stationary situation. This 
implies that F w (s; X) appears commonly in the stationary subclass in the KPZ universality. 
Our result above is a representation which does not include any resolvent kernel such as 
(1 — A)~ 1 (x, y) thus it is convenient for numerical calculation. In Fig. 2, we also depicted 
the picture for the long-time limit using (2.22). 

The derivation of the above theorem will be provided in Sec. 6. 

2.2 The two-point correlation functions 

The two-point function of the height is defined as 

C(x, t) = ((h(x, t) - (h(x, t))) 2 ) . (2.39) 
Considering the scaling property (2.4), we also introduce the scaling form 

*® = ^{*k)' C (2 - 40) 

where we put the overall factor 2 _4//3 and the factor 2 1 / 3 in the argument of C(x, t) following 
the convention in [45]. 

Moreover the two-point function of the slope of the height (d x h(0,0)d x h(x,t)) is also 
an important quantity since the height slope d x h(x,t) corresponds to the particle density 
in the interacting particle processes such as the ASEP. It has been pointed out in [45] that 
this is expressed as the second derivative of C(x, t), 

(d x h(0, 0)d x h(x, t)} = ^d 2 x C(x, t) (2.41) 

Thus we are also interested in the second derivative of the scaling function g' t '{y)- 

Noting that (2.39) is just the variance of h(x,t), we readily obtain the following result 
for g t (y): 

Corollary 4 

M = 2-/3 ( jf» ^l^^A is _ (|_~ t^f^jA (2 .42) 

where F Wit (s;X) is given in (2.22). 

Fig. 3 shows the function g"(y) for j t = 1,2,3. The dots in this figure correspond to 
the Monte Carlo simulation of density correlation in the WASEP 

(%,t?7o,o)wASEP - 1/4 (2.43) 

where r) x>t is the occupation variable, i.e. if the site x is occupied (empty) at time t, 
Vx,t — 1 (0) and ( )wasep means the averaging over both two-sided Bernoulli initial 
configuration and WASEP dynamics. 



11 




Figure 3: Stationary two-point correlation function g"(y) for the KPZ equation for jt = 
1, 2, 3. The dashed curve represents the one in the scaling limit g"(y). The dots indicates 
the Monte Carlo result of the WASEP, where parameters are the same as those in Fig. 2. 



3 Replica analysis 

3.1 Cole-Hopf solution, <5-Bose gas 

By the scaling of the KPZ equation, we find the height h Vt \ t D(x,t) of the KPZ equation 
with general parameters u, A, D and hi nix, t) with the specialized parameters v — ^, A = 
D — 1, are related as 

-^Kxd (x, = hi >ltl {a 2 x, aH), (3.1) 

where a is defined below (1.3), for details see [36]. Thus in the following we concentrate 
on the height hi 11 ix,t) and omit the indices |, 1, 1. 
Applying the Cole-Hopf transformation, 

Z{x,t) = e h{x ^\ (3.2) 

the KPZ equation (1.1) is linearized as 

dZ(x,t) ld 2 Zix,t) 



dt 2 dx 



2 



+ 7i(x,t)Z(x,t), (3.3) 



which is called the stochastic heat equation. Now we would like to remark on the regu- 
larization of the KPZ equation. The KPZ equation (1.1) is in fact not well-defined due to 
the interplay between the nonlinear term and the singular property of the noise in (1.1). 
The typical regularization is to interpret (3.3) as Ito type and to define the regularized 
KPZ equation by the inverse Cole-Hopf transformation of (3.3). This is called the Cole- 
Hopf solution of the KPZ equation [24]. For more detailed discussion, see [36]. Recently a 
new way of regularization without using the Cole-Hopf transformation (3.2) was proposed 
in 
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Eq. (3.3) can be interpreted as a problem of directed polymer in random media. Ap- 
plying the Feynman path integral representation and taking the average with respect to rj, 
we see the problem becomes that of the delta function Bose gas with attractive interaction, 



1 N a 2 1 N 

H ^— 2 Y.3^-- 2 Y, 5 ^-^ ( 3 - 4 ) 

3=1 j^k 

in terms of which the replica partition function (Z N (x,t)) can be written as 

(Z N (x,t)) = (x\e- H * rt \$). (3.5) 

Here (x\ represents the state with all N particles being at the position x and the |<I>) 
the initial state of the 5-function Bose gas. This is an important property of the replica 
approach to the KPZ equation [28, 29, 30]. The integrability of the Hamiltonian (3.4) 
allows us to calculate exactly the replica partition function. 

One can perform the average over the initial distribution given by the two-sided Brown- 
ian motion and the dependence of |<&) on x±, . . . , xn can be explicitly calculated as follows. 
First let us assume x\ < x 2 < . . . < xn- For x\ < . . . < xi < < xi + \ < . . . < x^, 1 < I < 
N, one finds 

(xi, ■ ■ ■ ,x N \$) 

= (exp + • • • + B-(- Xl ) + B + (x l+1 ) + ■■■ + B+(x N )]) 

l N-l 
= e v-T, l j=1 x j -v+j:f =l+1 x j JJ e -^(2J-2j+l)x, JJ e ±(2{N-l-k)+l)x l+k ^ ^3 

3=1 k=l 

where ( ) denotes the averaging over BMs B±(x) (2.1) and in the last equality we used 
the property of BM, (B(x)B(y)) = mm(x,y). Since we are considering a Boson system, 
we symmetrize the function with respect to xi, . . . , xn to have 

i 

(x ir -- ,x N \$) = e v -^=^ Xp ^- v +^=^ Xp ^Yl e ~^ i2j ~ 1)Xp(j) 
Pes N 3=1 

N-l N 

x Yl e ^-^)+D^ +fc ) JJ {xm _ xp{ . (3 . 7) 

k=l j=l 

Here Sn denotes the set of permutations of order N, Q(x) is the step function and we set 
xp(o) = 0. 

The eigenvalues and eigenfunctions of the 5-Bose gas can be constructed by using the 
Bethe ansatz [32, 33] (see also [30]). Let \^ z ) and E z be the eigenstate and its eigenvalue 
of H N , 

H N \V Z ) = E Z \V Z ). (3.8) 
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By the Bethe ansatz, the eigenstate is given as 
(x 1 , ■ ■ ■ , x N \t? z ) = C z ^ s & nP II ( Z PU) ~ z P(k) + - Xk)) ex P I * Zp (0^ 

PeSjv l<j<k<N \ Z=l 

(3.9) 

where C z is the normalization constant, for which a formula is given in (3.11) below. 

For the <5-Bose gas with attractive interaction, the quasimomenta Zj (1 < j < N), 
which label the state, are in general complex numbers. Zj (1 < j ' < N) are divided into M 
groups where 1 < M < N. The ath group consists of n a quasimomenta z'jS which share 

the common real part q a . Note that Y2^?=i n a — N. The quasimomenta in each group line 
up with regular intervals with unit length along the imaginary direction. Using q a and 
n a (1 < a < M), we represent Zj (1 < j < N) as 

a-l 

z j = Qa - ~ (n a + 1 - 2r a ) , for j = ^ n P + r «> (3-10) 

where 1 < a < M and 1 < r a < n a . The normalization constant C z in (3.9), which is 
taken to be a positive real number, and the eigenvalue E z are given by [30] 



a=l a l<j<k<N 1 J 1 / 

^ N M M 

o E z " = 9 E - oz E K - n «) • ( 3 - 12 ) 



2^ J 2 ^ a 24 

j'=l o=l a=l 



The completeness of these states has been proved in [53, 54, 55]. 
The replica partition function (Z N (x,t)) (3.5) can be written as 

/oo poo 
dyi--- dy N (x\e- HNt \ yi ,--- ,y N ){y u --- ,y N \<f>) (3.13) 
-oo J — oo 

Expanding the propagator (x\e~ HNt \yi, • • • , yjy) by the Bethe eigenstates of the 5-Bose gas 
(3.9), we have 



AT iV / poo M , oo > 

< ZA (*.')> £ ■rll/ <U / II ' EJ 

M=l j=l" / -° \" / -°°a=l n„=l/ 

x e-^I^K^Iy!,--- ,y N )(yi,--- ,Vn\^). (3.14) 

Now we want to perform the integrations over y^, (1 < j ' < N) in this equation, Us- 
ing (3.7) and noticing the symmetry of the eigenfunction (3.9), z \yp(i), ■ ■ ■ ,Up(n)) — 
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(®z\yi, ■ ■ ■ ,Vn), we find that RHS of (3.14) is represented as 

N AT\ N r M / „oo , oo \ 



x e- h - A {x\m z )(m z \ yi ,--- ,y N ) 

l N-l 

x e «-£i=iW- , '+£ J -= l +iW JJe-iPJ'- 1 )* JJ e |(2(JV-;-fc)+i)M +fc _ (3-^ 

i=i fe=i 

Here we find the necessity of the drift terms in (2.1). Note that due to the factors 
exp(— ^(2/ — 2j + l)yj) and exp(|(2(iV — I — k) + 1) we could not perform the integrations 
over t/k, k = 1, ■ ■ ■ , N prior to those over q a , a — 1, • • • , M if we did not introduce the 
drifts i>±. Now we perform the integrations of y 3 - (1 < j < N) by assuming that the 
drifts are positively large enough, which is the main reason why we introduced the drifts. 
Eq. (3.14) can now be expressed as 

N -. M / „oo , oo \ 

<*"<*.«» = E ./ It E * Eftl ^<*i».><».i»>«-*«. (3.i6) 

M=l a=l V^- 00 n a =l/ 

Here (a^z) is given by (3.9) with x 1 = ■ ■ ■ = x N = x and (&z\®) is computed as 

N ,00 

j=l J -°° 

= N\C Z {z* m - z* P{k) +i))jr f d Vl --- d Vl n e (- iz nm^-\^))y- 

7_n J Vl < — <Vl <0 1 



» JV— i 

/ dj/i+i • • • dy N ] [ e (-^(« +m )-"++3(2(^-'-"»)+i))« + « 

J0<Vl^T<---<VN 1 



PeSjv Z=0 J !'i<-<tt<0 m =i 

JV— i 

X 

'0<j/ i+ i<-<j/jv "Ji 

(-1)"JV!C7, £ sgnP J] {z* m -z* P(k) + i) 

PtS N l<j<k<N 
N I N-l 

X 



: j= iv- m+ i(-^po-) - v +) + m 72 

(3.17) 

Here we assume that Ejli(~^P(j) -m 2 /2 > and J2j ! =N-m+i(~'i z p(j) -v+)+m 2 /2 
for any m (1 < m < N). For example if we set v± > iV/2+max Q n a /2, the above conditions 
are satisfied. 
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3.2 Combinatorial identities 

For further analysis of the integrand in (3.16), we need two combinatorial identities for 
(x\^ z ) and {$ z \&). The first one is for (x\^ z ). One has 

sgnP J] (wp U) -w PW + if(j,k))=N\ J] ( Wj -w k ) (3.18) 

PeSjv l<j<k<N l<j<k<N 

for any complex variables Wj (1 < j < N) and f(j,k). This identity was derived as 
Lemma 1 in [37]. 

The next one for the term (fy z \Q) is 

Lemma 5 For any complex numbers Wj (1 < j < N) and a, 

^2 sgnP Yl ( W PU) - w P{k) + a) 
PeSjv i<j<k<N 

N l , N 

X 



J2( l) 1 f[ +v _}_ m 2 a , 2 II yN 

1=0 m=l ^J= 1V ^w; / / m= i 



j = JV-m+l 



— f + ) + m 2 a/2 



nLiK + ^ - IIi<j<fc<JvK- - wfc) ^ 3 1 

l\Z=i( w m + v- - a/2)(w m -v+ + a/2) 



Proof 

First we note that a similar identity, 

sgnP H (w P(j) - w P(k) + a) H — ^ = - — — — , 

Pes N i<j<k<N m= iE M - m+ i«'P(,)+™ 2 «/2 U m =i( w m + a/2) 

(3.20) 

which is the special case V- — > — oo, i>+ — >■ of the above (3.19), was already proved 
in [36]. By using this identity for both terms with v + and V- on LHS, one can rewrite LHS 
of (3.19) as 

TV 

5^(~ i ) n yi sgns n k-j _u, jv) n ( w k m -wk m ,) 

n=0 l<ji<...<j„<Af l<i<Z'<n l<m<m'<N-n 

l<fcl<...<fciV- n <7V 

{ji}7 =1 n{k m }^=<t> 

n N—n n . N—n . 



n n <«* - -~ + «) n n mfc ,.„ ++a/2 - ^ 



X 

1=1 m=l 



where S = {ji, ■ ■ ■ ,j n ,ki,--- , kjy^ n }. At this point we multiply both sides by rim=i(' u; »i + 
t>_ — a/2)(w m — v + + a/2) and then set v + — a/2, V- = 6 + a/2 without losing the generality. 
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Hence what we want to show is 



N 

$^(-i) n Yi s & ns n k-^v) n (vkm-vkj 

n=0 l<j 1 <...<j n <N l<l<l'<n l<m<m'<N-n 

l<hi<...<k N _ n <N 

{jj}P=in{*™}£I=r=* 

ra JV-n n JV— ra N—l 

x n n k - ^ + °) n w » n + ^ = n ( 6 - am ) n - o- 

1=1 m =l 1=1 m =l m=0 l<j<k<N 

Here we prove this by noting that both sides can be considered as a polynomial in a and 
b and arguing by induction that the coefficients (which are polynomials of {w±, . . . , wn}) 
for each b n a m in both sides are the same. First we observe that LHS is anti-symmetric with 
respect to {wi, . . . ,wjy}. Hence the coefficient b n a m is also anti-symmetric with respect 
to {wi, . . . , Wn}. Then all we should check is the numerical coefficients of terms Ylf=i w °j 
satisfying > • • • > C\ > 0. For instance when N = 2, the equality we want to show 
reads 



(w 1 - w 2 )(w 1 + b)(w 2 + b) - (w 1 - w 2 + a)w 1 (w 2 + b) + (w 2 — w 1 + a)(w 1 + b)w 2 

+ (wi — w 2 )wiw 2 = b(b — a)(wi — w 2 ). (3.23) 

We compare the coefficients of 1, a, b, ba, b 2 . For example for that for b, we see the coefficient 
of Wiw 2 is zero; for ab, we see the coefficient of w 2 is one. 

To consider general N, we introduce a graphical representation. Each graph corresponds 
to each term in LHS of (3.22). First we put n numbers . . . ,j n } to the left and N — n 
numbers {ki, . . . , k^-n} to the right. They are the vertices in a graph. Next if the term 
contains a factor from Wi — Wj we draw a line joining i and j (edge); we also put circle 
on % if the term contains Wi and j if Wj. In addition, if the term contains a factor coming 
from the last two products in LHS of (3.22), we encircle the number; hence we encircle all 
jVs. Finally by considering how each term can appear, we can put the sign to each graph. 
Because of the antisymmetry with respect to {w±, ■ ■ ■ ,wn}, we need to list up only the 
graphs with c N > ■ ■ ■ > C\ > 0. The example of these graphs for N = 2, 3 are given in 
Figs. 4 (a) and (b) respectively. By looking at these graphs, one can confirm that there 
are several graphs which give non-zero contributions and all others give zero. 

For showing the identity for general N, we use an induction. Suppose that we already 
showed the identity for iV by using the graphical representation. It consists of the several 
graphs corresponding to RHS of the identity (see e.g. the columns "6 3 ", "6 2 a" and "6a 2 " 
in Fig. 4 (b)), where the corresponding coefficient is nili^l 1 U P to sign, and all other 
graphs being grouped to vanish. Next we want to prove the identity for N + l. First notice 
that every graph which appears for iV + 1 is obtained by adding the vertex iV + 1 to some 
graph for N. If it is from a graph which can be grouped to vanish at N, then the term 
at iV + 1 can also be grouped to vanish. On the other hand, by adding iV + 1 to a graph 
corresponding to RHS of the identity at N, only four possible types of graphs at N + 1 
corresponding to the coefficient n^i 1 w j 3 with cn+i > cn(— N — 1) > cn~i(= N — 2) > 



17 



• • • > c 2 (= 1) > ci(= 0) can be produced. Fig. 4 (c) illustrates these four types. If we 
suppose the base graph at N represents the coefficient 11^=2 "^j" 1 °f b n a N ~ n , the types (i) 
and (ii) at N + 1 (produced by adding the vertex 11 N + 1" to the base one) correspond to 
the coefficient (— l) N+1 w N+1 Ylf=2 ^j" 1 an< ^ (^) N ' ujN+1 Y[f=2 w j~ 1 °f b n a N ~ n respectively, 
which cancel out each other. The types (iii) and (iv), on the other hand, correspond to the 
coefficient (-1) N+1 Uf^ 1 w]~ l o{b n a N+l ~ n and (-1)" oib n+1 a N ~ n respectively. 

Note also that in the type (iii), there are N choices in total for connecting the vertex N + 1 
and the other ones with lines. Thus we find that for m — 0, 1, • • • ,N 

C[b N+1 ~ m a m ] = (-l) N w% +1 C[b N - m a m ] + (-l^Nw^C^-^-^aT- 1 ], (3.24) 

where C[b n a m ] represents the coefficient of b n a m . This completes the proof. □ 

Using the identities (3.18) and (3.19) to (3.9) and (3.17) respectively, we get 

(x\* g ) = N\C z J] (zj-zje^i**, (3.25) 

l<j<k<N 

N _ j 

(»,!•) = mc, n W-on w + fa .^"_ fa+ + </2) - (3-aO 

l<j<k<N 1=1 v « I J\ I + / / 

Thus using these relations and (3.11), we find that the factor 1 3/^) | < ^ ) ) in (3.16) be- 
comes 

M 1 ^ Jv k--^-ii 2 fiH.r+.--i/2)(^+, + -i/2)' 

(3.27) 

We want to rewrite this equation in terms of g a and n a in (3.10). For the last factor of 
this equation, we easily find 

N 



11 (_^ r + w _ _ i/2)(^* + v+ - 1/2) 

. TT e j" Q g Q ^ TT * /g 23) 

II 11 (_ iqa + v _ + U na _2r))(iq a + v + + l(n a -2r)Y 



a=l r=l 



From (3.11) and (3.25), we know the remaining factors in (3.27) is represented by \ (0\^/ 
For this quantity, the following result was obtained in Appendix B in [30], 

f^a • ) TT | Zj ^A; | 



Koi^)i 2 =iV!n^ n 



a=l 71(31 l<j<k<N Zk l ' 



N\ tt \q a - Qp - j(n a - n p )\ 2 

11 \„ -n„-L<„ X„J|2- {d - 29) 
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Figure 4: Graphical representations for the coefficient of each term in LHS of (3.22). (a) 
and (b) represents the case N = 2 and N = 3 respectively. In (c), we depict the possible 
four types of graph at iV — 1 which comes from the graph at iV corresponding to RHS 
of (3.22). 
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From (3.27)-(3.29), we get 



t=\ a</3 
n a 

X 



n (-i Qa +v _ + ^ Ua - 2r))(iq a + v + + \{n a - 2r)) ' ^' 3 ^ 

We can further deform (3.30) to an expression in terms of a determinant by using the 
Cauchy's determinant formula, 

and a few basic properties of determinant. We find 

N / j \ M 



(x|^ 2 )(^|$) = 2 M iVlf[(t; + + ^-£) det f 



X 

q=1 r 



\7ij +n k + 2i(qj - q k ) J j k=1 
1 



I I e in a q a x 1 I 

H r A i Hfc. + u- + |K - 2r))(ig Q + v + + |(n a - 2r)) 

2 M iV! H(v + + v--e)Y[( dcoA det ( e -^K+n fc +^fe- ?fc )))^ =i 



£=1 a=l 
M n a , 

X JJ e*™"^ JJ — 



i (-«5a + v- + \{n a - 2r))(iq a + v + + |(n a - 2r)) 
= 2 M iV! JJ( U+ + u- - ^) II ( / det ^ W J> w *))ffc=i > 

£=1 a=l ^ ' 

(3.32) 

where in the last equality, we used a simple fact 



det(a 6j+6fc ) = det((a J a fc ) 6 0- (3.33) 



and the notation 



^injqjX-nj(ujj+uj k )-2iqj(ujj-uj k ) 



n 



\ (-iqj + v_ + \(rij - 2r))(?g j + v + + \{rij - 2r)) 



r + v _ + a ) r fa + v+ + ^) • [6 - 6 } 
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From (3.12), (3.16), and (3.32), we obtain an expression of (Z N (x,t)} in terms of the 
determinant, 



(Z N (x,t))e^ + — 

M / 00 r°° rln r°° \ 



N M M , N M / oo oo oo \ 

= E^nK+»-on e/ ^e-'^< 

M=l 1=1 a=l \n a = l J -°° ' Jo J 

fx 2 \ M 

x det ( e nj % K x (nj , qj \Uj,u h )\ 

V / j,k=l 

N N ) N M / oo oo \ 

= Ej^IK W + + W -- € )II / ^ E ^n,iV 
M=l ' £=1 a=l V 70 ra a =l/ 

x det f /°° ^ e "^-t"^ 2 +^"f K x {nj, q; u,, u k )\ . (3.35) 

Here we would like to discuss the condition on the drifts v±. In the above equation, we 
notice that the condition v± > N/2 + m&x a n a /2 discussed below (3.17) ensures that all the 
poles q = —i(v- + (rij/2—r) and q = i(v + — (rij/2—r) with r = 0, 1, • • • , rij in the integrand 

e nj 2t~2 njq2+ 24 n iK x (nj, q; ujj, Uk) lie on the lower and upper half plane respectively. Thus 
we can relax the condition to v± > max a n a /2, where the above property of the poles is 
retained. 

4 Generating function and height distribution 

4.1 Generating function 

Now we define the generating function G 7t (s;X) as 

G lt (s;X) = J2 ^~~~Jy~\ ( ZN [HX,t))e N ^ x2 = {e^ HtW ~'\ (4.1) 

N=0 

where in the last equality, we used (2.4). Substituting (3.35) into this equation, we get 

G lt (s;X) 

oo N N (_- ltS \N M ( poo °° \ 

= i + e riK + «- - o e t ^ n / ^ e 



Ml HI /„ 

JV=1£=1 M=l a=l V u n a =l/ 



x det ( r ^ e ^^»^+£»?/^(n,-, g; Uj , u k )) " . (4.2) 

\J-oo K J j t k=l 
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4.2 Shifting 

In the previous studies for the narrow-wedge [30, 28] and half Brownian motion initial 
condition [36], the generating functions corresponding to (4.2) were already described by 
the Fredholm determinant, which is an important achievement in the replica analyses. 
Unfortunately, however, this is not the case in (4.2) due to the factor [T^LiC 1 ^ + v - ~ j)- 
In order to resolve this problem, we introduced h(x,t) in (2.6). Here we discuss some 
general properties of two random variables h and h related as 

h = h + X , (4.3) 

where x is & lso a random variable independent of h. Note that at this stage, h and h are 
just two random variables and do not have a specific meaning such as height. In addition, 
right now we do not assume a specific distribution for x- Let us define 

Z = e\ M(0 = £ = <e^> , F(s) = F[h < s]. (4.4) 

N=0 

and also Z, M(£), F(s) in the same way as (4.4) except that h is replaced by h. We find 
that the following relation holds 

{zN) = w~f N) (45) 

M(() = (4-6) 
and noticing the fact = d£M(i£)e~^ s , we also see 

F M - (") 

where g(y) = (e yx ). 

In (2.6), we chose x as the random variable such that e x obey the inverse gamma 
distribution with parameter v + + t>_. The probability density function of the inverse 
gamma random variable is 

Pe(x) = ^-^e-V^o, (4.8) 
where 9 > is a parameter and the function g(£) in this case is given by 

9(0 = j™xtp e {x)dx = -L J\^ e -^dx = (4.9) 

In particular, we notice that if we set 9 — v + + v- as in (2.4), g(N) becomes i/YljLi( v + + 
v- - j)- 
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4.3 Generating function for the generalized height 

We define Z(x,t) and G lt (s;X) as (3.2) and (4.1) respectively with h(x,t) replaced by 
h{x,t). From the relation (4.5) and (4.9) with 9 = v + + t>_, we see that (Z N (x,t)) is 
related to (Z N (x,t)} as 



1 

ni=i( v + + v - ~ j) 



(Z N (x,t)) = (Z N (x,t))- ¥r - / ■ -. (4.10) 



Hence G lt (s;X) is calculated as 

°° (—p-lts\N 

G lt (s;X) ■= £ m (Z N (^X,t)) 

N=0 

co N , ~ ttS \N M ( roo oo \ 

JV=1M=1 ' a=l V u n a =l/ 



x det QH ^e n ^ x2 -^ 2+ ^K Hx {n ])q -u ^^ 



M 



j,k=l 

(4.11) 



Note that this generating function does not include the term n^i( l '+ + V- — £) in (4.2). 
In this case, we find the Fredholm determinant expression, 



°" (_-i\M 1V1 poo 
M=0 ' a=l J ° 

x det (jrt-iy 1 - 1 [°° ^e"^* x2 ^ 3 "^ 2+ ^"?^" s K 27t2x (n,g;^,o; fc )) . (4.12) 

\n=l J ~°° 71 J j, k =l 

Shifting the variable q to q + iX/^t, we get 

00 f-l^ M poo M 

G lt {s-X) = J2 ~lJrU / ^ det , (4-13) 

M=0 ' fc=l ^° J ' 

where the kernel K VjW (ojj,Uk) is written as 

oo /*oo J 3 

K vw (u j: u k ) = V(-l)"- 1 / ^ e -™K+^)-2^K-^)- 7 3„ 9 2+2|„3- 7tns 

x r(-ig + w -z)r(iq + v-*) 

r + iw + f ) r (ig + v + § ) ' 

with t>,w > n/2. In (4.13), we keep the contour of q on the real axis by changing the 
condition for the drift to v + — X/^ t > n/2, v _ + X/^ t > n/2. 
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In the following, we provide another representation with the kernel using the deformed 
Airy functions Aip(a,6, c, d) defined in (2.14). We use the following two relations, 

Lemma 6 

(a) We set q el and m, n > 0. When a,b > n/2, we have 



r (-iq + a- f)T (iq + b- i) m^i _ f°°^ iX , rr f ni 1 



r(-iq + a+2)r(iq + b+%) 



J dyM™ (y, -iq + a,iq + b^j e mny , 



(4.15) 



where 



Ai rr (abed)- — f dzc iaz+iZ " /3 T (' bz + c ) T ^ bz + d ) u i 6 ) 

Ai rr [a, b, c, d) - ^ ^ dze p( _ + c)r( _ + rf) (4.1b) 



with T p defined in (2.17) below. 

(b) For a,b,x 6t and w > 0, we /iawe 

1 Z" 30 

— / rfpAipp (p 2 + f , w, —iwp + a, iwp + b) e ipx 

= ^r Ai r (2~^(v + x),2^w,a,b^ Ai^ ^2~i(u - x), 2^w, b, a) , (4.17) 

where Ai r (a, b, c) and Air(a,b, c) are defined by (2.18) and (2.17) respectively. 

As will be discussed in Sec. 7, these relations (4.15) and (4.17) play a crucial role in 
the replica analysis. These are two-parameters (a, b) generalization of the ones used in the 
previous studies of half-Brownian motion [36] and the narrow wedge [28, 30] initial data: 
As either a or b goes to infinity in (4.15) and (4.17), they reproduce Lemma 6 (a) and (b) 
in [36] respectively. On the other hand, in the case where both a and b go to infinity, they 
converge to the ones appeared in [28, 30]. 

Proof 

(a) RHS of (4.15) is written as 



.3 



1 f 00 f T (i^- -iq + a)T (i^- +iq + b) ., . , 

— dy dz / 2m ( 2m ^ j r e^^T. (4.18) 

2tt J-oc Jn+imn T {-i^ - iq + a) T {-i^ + iq + b) 

Here we used the fact that the contour of z can be deformed from to R + imn since the 
imaginary parts of the poles of the integrand 2m{ia — q + ir) and 2m{ib — q + ir), where 
r = 0, 1,2, ••• are larger than mn when the condition a,b > n/2 is satisfied. In this 
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equation, we change the variable z on RHS to y 2 = z — imn and get 



2k J -00 J-00 rH g + a + f-^)r( J9 + Hf-^) 



^T(-^ q + a + ^-^M ^ )T(^ q + b+■l-^^) 
r (-iq + a- f)T{iq + b- f) 



/oo 
-00 



(4 19) 

r (-iq + a + f)T (iq + b + f) ' V ' 7 

(b) LHS of (4.17) reads 

1 r . / dzei(p i +v)z+iS i +ipx r (n^z - g) + q) r (n^z + p) + 6) 

(2ir) 2 J_ OB ab r(-iw(z+p) + a)r(-iw(z-p) + by 



(4.20) 



By applying the change of variables p — {z\ — z 2 )/2 2 ^ 3 and z — (z± + z 2 )/2 2 ^ 3 , we obtain 
the desired expression. □ 

Applying (a) of the lemma with a = v-+X/ r y u b = v + — X/^ t and m = 7j/2 2 / 3 to (4.14), 
one has 



v+ it 1 ' it 



^ J-OO * J-OO V 23 7t 



+ £ >i5 + u+ _^) e -2«/(^-^) 



It It 



. 2 



x e -n{.Uj+Uk+ltQ +1ts~2 3-ytyi) (4 21) 

So far we have imposed the condition for the drifts v+ — X, v- + X > n/2 (see (4.13) 
and (4.14)). At this point, since RHS of (4.15) is well-defined for all o, & £ 1, we can 
relax it to the region such that the generalized height (2.6) under consideration is well- 
defined, i.e. v + + v _ > 0. Changing the variables q, coj and yi to p — 2 1 / 3 7 t g, £j = 
2j~ 1 uj j and y = yi/2 2 / 3 -p 2 /2 2 ^ - + £ k )/2, we see K^x/^v.+x/yt^j,^)^ = 
Kv+-x/ft,v-+x/*ft(tj> £k)d£j, where 

2 ^ Z" 00 ^ Z" 00 ^ A-rr/ 2 , o- , + & 1 -ip + 2i-f t w ip + 2i~f t v 



2§ z" 00 z" 00 / 2 

tk) = — J dy J dpMll [p 2 + 2ay + 



2s 2s7i 237 t 237 t 



x e 2 i ^(_i)"-i e -w»(«-w) 

n=l 

^Atf (fc + y, - , Ai? + y, -,w,v) elty + elts - (4.22) 

where in the second equality we applied (2) of Lemma 6 to this equation with v = 2 2 l 3 y + 
(0 + 6)/2 1/3 , w = l/(2 1 /3 7t ), a = v _ + X / 7t , & = X/ 7t , * = -(fc - 6)/2 1/3 . Note 
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that although the sum Y^=i(~ l)™ _1 e _7 *"( s_2/ ) in the first equality become divergent for 
s — y < 0, the last expression e lty /{e ltV + e ltS ) in the second equality is well-defined in 
arbitrary s — y G R. As we will comment in Sec. 7, this divergence originates from the 
ill-defined nature of the generating function and this kind of analytic continuation has 
been necessary for the previous studies on the replica analyses [28, 30]. Using (4.22), we 
eventually obtain 

G 7t (s;X) = det (l - P K x >v + *P ) ■ (4.23) 

v ^ -ft ' -ft / 

4.4 Proofs of Proposition 1 

Eq. (2.9) follows immediately from (4.7) and (4.9). For the derivation of (2.10), we use 
a formula given in [28, 37, 36] which transform the generating function G lt (s;X) into 
the height distribution F v+iV _ it (s; X) (2.5). By using the Fredholm determinant (4.23), 
F 7t (s;X) can be expressed as follows [28, 37]. 

/oo 
due- eltis - u) g Jt (u;X). (4.24) 
-oo 

Here 

= ^ (det(l - PoK + _± >v _ + ± P ) - det(l - PoK; + _^ v _ + ±P )) (4.25) 

where _x_ x_ ( x , y) is the kernel (4.22) in which the term e~ s is replaced by — e u ± it 

v+ a ' v ~~^ -ft 
with e > being infinitesimal. 

Using (4.22), and the relation l/(x±ie) = V(l/x) =f iir5(x), where V denotes the 

Cauchy principal value, we can easily find in (4.25) is represented as 



/OO J 



1 


X 


X 


— ,v+ - 




,^- + — 


It 


It 


It 


1 


X 


X' 


— ,v+ - 




V- + — 


It 


It' 


It, 



1 X X 

■ + u i i v - H ,v+ 

It It It 

(4.26) 

Substituting this expression to (4.25) and using basic properties of determinant, we even- 
tually arrive at the expression (2.11) in Proposition 1. 



5 Stationary limit 

In this section, we provide the proof of Theorem 2, which represents the exact height 
distribution in the stationary limit (2.19). This is obtained by rewriting the representation 
in Proposition 1 in a more suitable form for taking this limit. 
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We first focus on the part det (l — P u (BF /t — P£j)P u ) in (2.11). Introducing the scaled 
drifts ou± as v± = cu±/7t, one finds that it is written as 

J-oo V It It It J \ It It It J 

(5.1) 

where Cf\y) is defined after (2.27). Thus it is expressed by a product of integral kernels 
P u {Bl - PldPuiii, 6) = X 1 X 2 (^ 1 , 6) where 

v) = P«(6)AiF + V, -, , (5-2) 

V it it it J 

X 2 (y, 6) = ($%)A$ U 2 + y, -, <±±^, ^^-) P tt (6) (5.3) 

V 7t 7t 7t / 

By using the relation of the Fredholm determinant det(l — X X X 2 ) = det(l — X 2 Xi), we 
have 

det (1 - P u (Bl - P^P U ) = det(l - X 2 X 1 ) (5.4) 
where the kernel X 2 X\ is expressed as 
X 2 X!(y 1: y 2 ) 

= C t w (j/i) / d£Ai£ [yi + £, -, , — Aif; y 2 + - , 



7t 7t 7* / V It It It 

r {S) (ll \ r r r ( "+~ x ^ +{)r f ^-+ x +^ + 1 ■ 
= ^*i§2 / tfal / dx2e iMyi+u )+ iMy*+u) +U 4+4) V ^ ^ ^ £ 

( 27r ) Jr i(w+ _ x) Jr i(w _ +x) r + 1) r ( "+~*~" 2 + lj 

—1 w_ + X — ix 1 u + — X — ix 2 

x — . (5.5) 

i(x\ + x 2 ) cj + — X + ixi U- + X + ix 2 

Noticing the relation 

1 b-ixx a-ix 2 1 , . 1 1 , . 

-(a + b)- (5.6) 



i{x\ + x 2 ) a + %x\ b + ix 2 i(x\ + x 2 ) %X\ + a ix 2 + b 

we can further rewrite it as 
X 2 X 1 (y 1 ,y 2 ) 

= cf\ yi ) T^Ai? Le-.l + ^H^ 

A V 7t 7t 7t 

x Aij; h/2 + e, -, 1 + — — , 1 + 

V it it it ) 

(yi)B LU _ + x,uj + -x,u{y2), (5.7) 
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where B a ^ u {x) is represented as 

B ahu ( x ) = -^ dz— >^ f (5.8) 

2tt 7 ria tz + a r ^ + i ) 

and we find this can be rewritten as (2.20). 
Thus we have 

det (1 - P u (Bl - P^Pu) 

= det (l - A^ + _ x ^_ +x ^ (l - (u+ + w_) J dxipA^C^B^x^+x^^B^+x^x^x) 
= (o; + +o;_)det (l-A^_ x ^_ +x ) 
x ( - — — — — - / dx(Cf ] 'B UJ+ _ Xt0J _ +XtU )(x)B IJJ _ +Xt0J+ _ XjU (x) 

V U + + U - J-oo 

+1- / dx(AWp Alt) C®B u+ _ x „_ +XtU )(x)B u -+X,u) + -X,u{x) ) 

J — oo / 

= K + (det (l - 4»-x*- + x) C-x,.. + x(u) ~ det(l - " • 

(5.9) 

Here and are defined in Theorem 3 and p A (6)(x,y) := (1 — A^)~ 1 (x, y). We 

can also rewrite the other term, det (l — P u B^ t P u ) in (2.11) in the completely parallel way. 
We have 

det(l-P u £?;P u ) 

= (oj + + u-) (det (l - A^-x^+x) L u+ - X ^+x(u) - det(l - A; + -x,u;_+x - A^+-x,u;_+x)) 

(5.10) 

Using these forms (5.9) and (5.10), we finally obtain the desired expression (2.22), 
F V)tt (s,X)= lim F^+^ t (s;X) 



I'f^- • ') / d/ds . , 

hm — ^ - 1 + i -FV^isjX) 

p / , ^ \ y UJ + + UJ^ J -ft ' -ft ' 

\ J 



p / UJ++UJ- _|_ J 



d/ds r _ e 7t(s-u) 



lim — ^ I .(s) / <iwe e Q^fw) 

V=- + r (" ++ ^ t +d/ds + i) V ~ ~* + y_ 00 

rn ^ . , f d^^K.^X) - (5.11) 

r(i + 7 t rf/rfs) 
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Here u Wjt (u;X) and is^ t (u;X) are given by (2.22). In the last equality, we used the fact 
lim w _^+w_ Ff± f__ t (s; X) = 0, which follows from the divergence of the term x m (2-6) in 

UI+=W -ft ' "ft ' 

the stationary limit u + + w_ — > 0. 

Finally we derive the representation in (2.28). By using B^ bu and B^\ u defined 
in (2.29), L ab {u) is written as 

1 „ l. poo 

LaM = -^y- - J JxC t (x)BU u (x)BU u (x) 

/CO POO 
dxC t {x) (b^ u {x)bV u {x) + B«»BM - / d *C t (x)B^)Bl 2 Ux). 
-oo 7 J —oo 

(5.12) 

Thus it is enough to show that in the limit b — > —a, the first two terms of this equation 
converge to —27/74 + w — a 2 — 1 in (2.28) since the convergences for the remaining terms 
are obvious. Using the formula 

/ C t (x)e- yx dx = - cot ( — J (5.13) 

J -00 it Kit J 

for < Re(y) < 7 t , we can further rewrite the term j^dxC t (x)B^l u (x)B^l u (x) as 

- r dxC t (x)B^ u (x)B^ u (x) = - e ~^~^ cot (£±V) (5.14) 

J -co r(^ + ij 7t v 7t j 

for < Re (a + b) < 74. Considering the following properties, 

r(l +x) = 1 --fx + 0(x 2 ), cot(x) = --^ + 0(x 2 ), (5.15) 
for x — > with 7 being Euler's constant, we easily find 

lim l ~ a ~u h ~ r dxC t{x)B { al u {x)B£l u {x) = - 2 -l +u -a 2 -l. (5.16) 

b^-a a + b J_ 00 ' ' ' ' It 

We notice that although the function L a ^ b (u) should be valid for < a + b corresponding 
to the condition for the drifts v± discussed below (4.21), the more strict condition < 
a + b < 7< is necessary for (5.12) in order to apply the relation (5.13) to (5.14) and other 
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terms in (5.14), 



r azc, w <>„ (,)<„(,) = ^- / * cot (az",) trrr r+1 l 



g— a!i+ — 



x — r - ^ (5.17) 



/ d^l / 


dz 2 COt I 


/R+ici •/] 


R+ic 2 ^ 



(2vr) 2 7i J M+ici J R+iC2 V 7t 

el(zi+Z2)u+t( , l 3 +2 3 ) r («±hi + 1) r (*±ja + i) 

(a + izi)(& + iz 2 ) p ^ h=m +ljr ( ^-' iZ2 + ij ' 

where c (ci, c 2 ) satisfies the conditions b < c < b + j t and < a + c < j t (a < c\ < a + j t , 
b < c 2 < b + 74 and < ci + c 2 < 7t). An advantage of this expression is that we can take 
the limit b — > —a for L a ^{u) as in (5.16). Moreover RHSs of (5.14), (5.17) and (5.18) are 
analytic also for a+b < and thus so is L at b(u) (5.12). Noting the relations (2.9)-(2.11) and 
considering RHS of (2.11) is represented by using (5.9) and (5.10), we find F v+:V _ :t (s; X) 
is analytic also for the case v + + v _ < 0. 

In the case j t < a + b, on the other hand, another expression for B a ^ u (x) 

f°° ( 1 b a\ 

B aA u(x)= d\e- aX Ai r r lx + u-\,-,l + -A + - , (5.19) 
Jo V It It It/ 

is convenient. Using this, L a;b (u) is calculated as 

l-a-b vr f f f-i( Zl + Z2 )^ e i^hH(4+4) 

La,b{ U ) = ~T ,r, \o I dz l / dz 2 COt' 



a + b (27r) 2 7 t J R+igl J R+i52 \ *y t J (a + iz 1 )(b + iz 2 ) 

Y I a+iz 1 _)_]_] p [ b+iz2 i ^ \ 

r(^ + i)r(^ + i) 

where c\ and c 2 satisfy the conditions ci < a, c 2 < b and < c.\ + c 2 < It- 

6 Long-time limit 

Here we give a proof of Theorem 3. Noticing the limiting behavior 

lim e~ eltx = l x < , (6.1) 

t— >oo 
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we find the long-time limit of Theorem 2 becomes 

F w (s;X) = - / du(u w (u;X)-u^(u;X)) = u^(s;X)-u w (s;X), (6.2) 

J S 

where v w (s;X) is defined by (2.32) and Uw\s;X) is given by 

X) = det (1 - A<») C^_ w+X (s) + det (l - A® - V^_ w+X ) . (6.3) 

Here A {s) and V[ '_ n are defined in the same manner as A (2.33) and T> a -a (2.35) respec- 
tively with Pq(x) replaced by Po(x) — 5(x). The function C^_ a (s) is defined in terms of 

C a ^ a (s) (2.34) and B a , s (x) (2.36) as C®_ a (s) = £ ,_ (s) - B a , s (0)£_ a , s (0). Thus for the 
proof of Theorem 3, it is sufficient to show 

v { J\s;X) - v w {s-X) = ^ w {s\X). (6.4) 

For this purpose, we focus on the expression in Proposition 1 and show that (2.9) has 
two equivalent forms corresponding to LHS and RHS of (6.4) in the simultaneous stationary 
and long-time limit (see (6.16) and (6.20) below). First, taking the long-time limit in (2.9) 
with v± = oj±/ r y t , we have 

F u+ p_(s;X) : = lim F w+ / 7t , w _/ 7t ,t(s; X) = ( 1 + j F w+iW _ (s; X), 

t— s-oo y cj + + cj_ as J 

POO 

F w+>w _(s;X) := lim F UJ+huUJ _ hut (s\X) = 1 - / dug(u;X), (6.5) 

where g(w;X) is defined as the long time limit of g lt (u;X) (2.11) with v± = oj±/^ t and 
has the expression as follows, 

g(s; X) = det (l - P s (P+_ - P/f )P S ) - det (1 - P S B + _P S ) . (6.6) 

Here B> + _ = lim^^ P^ and = lim^^ P£ are given by 

POO 

P+-(6,6)=/ dyAi + _(ei + y,o; + -X,a;_+X)Ai + _(ei + y,w_ + X,a; + -X), (6.7) 
Jo 

= Ai+_(ei,o; + +X)Ai+_(6,w- + -X), (6.8) 

where 

Ai + _(a,6,c) = lim Ai^ ( a,-,-,-} = [ dze iza+i ^^^. (6.9) 

V 7t 7t It J Jr ic c + iz 

Next taking the stationary limit {oj_ — >■ — u + , u + = w) of F LU+tUJ _(s; X) (6.5), we have 

lim F u+jU _(s;X) = lim g(s;X), (6.10) 
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where we used the fact lim w _^-u; + F u Ul _ (x, X) = as in (5.11). 

u>+=w ' 

Now we notice that RHS of (6.10) can be rewritten in two ways. First, we use the 
results (5.9) and (5.10) in the previous section for the expression of g(s;X) (6.6) . By 
taking the long-time limit for (5.9) and (5.10), we see 

det (1 - P U (B + ^ - P^-)Pu) = (w+ + (a; X), (6.11) 

det (1 - P U B + _P U ) = (u+ + u-)u u+ ^_(s;X), (6.12) 

where 

v a , b (s; X) = det (1 - A) C a+x , b -x(s) + det (1 - A - V a+X ^ x ) , (6.13) 

A(x, y), C a> b(s), V a ^(x, y) and B a ^ s (x) are defined by (2.33)-(2.36). The function i/i+ )W _(s; X) 
in (6.11) is defined in the same way as (6.13) with P replaced by Po(x) — S(x). Note that 
£ a ,b(s) has the following expression 

1 _ " 3 + >3 -(a+b)s poo 
£a,b(s) = — b 1 + jf dx [BQ(x)B^(x) + BU(x)B$(x)) 

poo 

- / dxBQ(x)I$(x), (6.14) 
Jo 

which corresponds to (5.12) and (5.14). From (2.33)-(2.36), one easily finds 

hm u u+ ^_(s;X)=u w (s;X), lim (s; X) = u^(s; X). (6.15) 



W+=W W-\-=W 



Thus from (6.6), (6.12) and (6.11), we find that RHS of (6.10) is rewritten as 



lim g(s;X) = u < £\s-X)-u w (s-X). (6.16) 



-w+ UJ , + UJ 



For the second expression for RHS of (6.10), we apply to (6.6) the relations 
det (1 - P s (B + _ - P£r)P s ) 
= det (1 - P S B + _P S ) (l + J dx ((1 - P s B+_P s )- x P s Pirp^ ( x , x fj , (6.17) 

/oo 
dx ((1 - P a B + .P 8 )- 1 P 8 P+-p 9 ) ( x ,x) 
■oo 

= -^trlog(l - P S B + -P S ) = -^-logdet(l - P S B + _P S ). (6.18) 
as as 



As a result we find g(s; X) (6.6) is expressed as 

d 

ds v + v v ' ds 



g(s;X) = -det(l-P s B + _P s ) = (u+ + u.)—u u+tU _(s;X) (6.19) 
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where in the second equality, we used (6.12). Hence using (6.15) we obtain the second 
expression for RHS of (6.10), 

lim 1 g{s-X) = ^-v w {s-X). (6.20) 

u>+=w 

From (6.16) and (6.20), we find that (6.4) holds. 

At the end of this section, we briefly comment the equivalence between our representa- 
tion in Theorem 3 and the one obtained in the studies of the PNG model [7]. From (6.5) 
and the first equality in (6.19), we have the following expression 

F„ +tUJ _(s;X) = lim F w+huUJ _ ha (s- X) =(l + ± det (1 - P S B + _P S ) . (6.21) 

Noticing that B + _(x) (6.7) can be rewritten as 

poo 

B + _(x,y) = / dXM(x + X)M(y + X) + (u + + u-)B u++x ,s(x)B u .-x,s(yl ( 6 - 22 ) 
Jo 

we find that this expression is equivalent to the one in Theorem 5.1 in [56], which is the 
limiting height distribution function in the PNG model with external sources. The relation 
between (6.21) and the expression in [7] was also discussed in Proposition 5.2 in [56]. This 
ensures that our result in Theorem 3 is equivalent to the representation in [7] in terms of 
the solution to the Painleve II equation. 



7 Conclusion 

In this paper we have considered the one-dimensional KPZ equation (1.1) in the stationary 
situation, which corresponds to the growth from the two-sided Brownian motion initial 
condition (1.3). Combining the Bethe ansatz results of the one-dimensional attractive 8- 
Bose gas with some techniques developed in this paper, such as the combinatorial identity 
and the shifting procedure, we have obtained the compact representation of the probability 
distribution of the height. Moreover the variance of this distribution provides the exact 
solution of the two-point correlation function. This function is a most fundamental quantity 
in statistical physics characterizing the space-time correlation and even if we focus only on 
the KPZ equation, it has been discussed by various approaches [57, 58, 59, 60]. Our exact 
solution in this paper should help us understand the nontrivial space-time correlation for 
the KPZ equation and universality in a deeper way. 

An important point for getting these exact results is to find out the Fredholm de- 
terminant structure for the generating function (4.23). However as pointed out in the 
introduction, the generating function is in fact a divergent sum. This can be confirmed by 
noticing from (3.12) that the ground state energy of the A^-particle (^-function attractive 
Bose gas is proportional to — A^ 3 and thus the N replica partition function (Z N (x,t)) (3.5) 
is proportional to e^ 3 . We avoid this difficulty in the following way: Using Lemma 6 
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(a), we rewrite e N in terms of e N and then we take a partial sum with respect to N 
in (4.22). (Note that some kind of analytic continuation was used in this equation.) Al- 
though the method involves some tricky procedures, we can say that our final results are 
fairly persuasive because of the following reasons. 

1. The rigorous version of the replica analysis has been recently studied in a certain 
interacting particle process called the g-TASEP of which an appropriate continuum 
limit provides the KPZ equation [39, 42]. 

2. As explained in Sec. 2, our result in Proposition 1 includes the exact solutions with 
narrow wedge [18, 19, 20, 21, 28, 30, 31] and half BM [25, 36] initial data as special 
cases. For both cases, the rigorous analyses by using the ASEP [18, 19, 20, 21, 25] 
has been done and the equivalence between these results and the ones by the replica 
method [28, 30, 31, 36] has been confirmed. 

3. In LHS of (4.17), the new gamma function factor appears in addition to the factor 
e m n / 3 , which is the origin of the divergence of the generating function in the case 
of narrow wedge initial data. However from the asymptotic property of the gamma 
function, we easily find the gamma function factor does not cause a singular behavior 
especially for large n. Thus the divergence of the generating function comes only 
from the factor e m n / 3 , which is the same situation as the narrow wedge initial data. 
Moreover as explained in (4.17), Lemma 6(a), which plays a crucial role in the replica 
analysis, is a natural extension of the cases for the narrow wedge and half BM initial 
data. 

To understand the regularization procedure in the replica analysis more deeply is an in- 
teresting future problem. We expect that the findings in the replica approach to the KPZ 
equation should lead to the promotion of understanding the replica method in general 
disordered system. 
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